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^ ■ Abstract 

We determine 0(a) boundary improvement coefficients up to 1-loop level for the Schrodinger 
Functional coupling with improved gauge actions including plaquette and rectangle loops. 
| These coefficients are required to implement 1-loop O(a) improvement in full QCD simula- 

. tions for the coupling with the improved gauge actions. To this order, lattice artifacts of step 

"j—? | scaling function for each improved gauge action are also investigated. In addition, passing 

through the SF scheme, we estimate the ratio of A-parameters between the improved gauge 
I ' actions and the plaquette action more accurately. 

1 Introduction 
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The MS-scheme now becomes the standard renormalization scheme for the definition of the 
strong coupling constant. The measured coupling constant a s in some experiment at relatively 
high energy is converted to ctj^g at some representative scale by perturbation theory. The 
current world average of such estimates gives a-^{mz = 91.19GeV) « 0.11. Lattice QCD 
calculations, on the other hand, have a potential ability to determine the the strong coupling 
constant from the experimental inputs at low (hadronic) energy scale. In order to compare 
the coupling constant obtained at low energy by the lattice calculations with ot-Tm obtained at 
high energy, the Schrodinger Functional scheme has been proposed by the ALPHA collaboration 
PQ, and the scheme is shown to be successful. At present, the results on the running coupling 
constant of two massless flavor QCD are reported jHIS]. 

In the real world there are three light quarks. QCD simulations including three dynamical 
quark effects are thus required to understand the low energy QCD dynamics. Our ultimate goal 
is to estimate from Nf = 3 QCD simulations. Recently, CP-PACS/JLQCD collaborations 
have started the iVf = 3 QCD simulation employing an exact fermion algorithm developed for 
odd number of quark flavors 00 IE]- In particular, noteworthy results have been obtained in 
[3]: There exists strong lattice artifacts associated with the phase transition in the Af = 3 QCD 
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simulation with the combination of the plaquette gauge action and 0(a) improved Wilson quark 
action, while such lattice artifacts are absent for the renormalization group (RG) improved gauge 
action. Hence, the collaborations have decided to adopt the combination of the RG improved 
gauge action and 0(a) improved Wilson quark action for the Nf = 3 QCD simulations to obtain 

a MS- 

As a first step of our program, we study a Schrodinger Functional coupling with improved 
gauge actions in perturbation theory. In particular we perturbatively calculate the 0(a) bound- 
ary improvement coefficients at 1-loop level with improved gauge actions for the pure SU(3) 
gauge theory. Combining this with the 1-loop results for the fermionic sector we determine 
the 0(a) boundary improvement coefficients at 1-loop, which can be used for the dynamical 
quark simulations in the future. 

The rest of this paper is organized as follows. In Sect 2, after a brief reminder of the 
Schrodinger Functional scheme and its extension to improved gauge actions, we specify an 
action used for latter calculations, and discuss the 0(a) boundary counterterm. In Sect 3 the 
Schrodinger Functional coupling constant is defined and formula for a determination of the 0(a) 
boundary improvement coefficients are given. The 1-loop computation is outlined in Sect 4, 
and the results of the 0(a) boundary improvement coefficients are summarized in Sect 5. Our 
conclusion is given in the last section, together with a discussion on the lattice artifact of the 
step scaling function. 

2 Preliminaries 

2.1 Schrodinger Functional 

It has been shown by ALPHA collaboration that the Schrodinger Functional (SF) scheme is a 
powerful tool to probe energy evolutions of some physical quantities and to compute improve- 
ment coefficients as well as renormalization constants. In the SF scheme, the theory is defined 
on a finite box of size L 3 xT with periodic boundary conditions in the spatial directions and 
Dirichlet boundary conditions in the time direction. In the pure SU(3) gauge theory with Wilson 
plaquette action S[U], the partition function in the SF scheme (in the case that T = L) is given 
by 

Z = J D[U]e- s[u \ (2.1) 

where link variables U{fx, x) for gauge fields satisfy boundary conditions 

U(x, k)\ Xo=0 = exp{aC}, U(x, k)\ XQ=L = exp{aC"}. (2.2) 

Here a is lattice spacing, and C, C are diagonal traceless matrices, which depend on background 
field parameters r/ and v [Sj. It is shown [T| that the minimum of S[U] is given by the lattice 
background field U(x,fi) = V(x,fi), where 

V{x,0) = l, V(x,k) = V(x ), (2.3) 



2 



with 



V(x ) = exp{a6(x )}, (2.4) 
b(x ) = ~[(L-x )C + x C]. (2.5) 

This background field represents a constant electric field. 

An extension of the SF scheme to improved gauge actions was first considered by Klassen 
[2]. The transfer matrix construction ^Qj was adopted in the discussion. In this formulation, 
Each boundary consists of two time slices, to achieved the tree-level 0(a 2 ) improvement. 

In this paper, however, we adopt the formulation proposed by Aoki, Frezzotti and Weisz 
jllj . where each boundary consists of only one time slice and the tree- level 0(a) improvement 
is achieved: Dynamical variables to be integrated over are independent on the form of the 
action, plaquette or improved ones, and are given by the spatial link variables U(k, x) with 
xq = a, • • • , L — a and temporal link variables U(0, x) with xo = 0, • • • , L — a on a cylinder with 
volume L 3 x L. This formulation is implemented more easily in numerical simulations. 

The background field in ea. (|2.3|) gives the local minimum of the theory in both cases |5| and 
jllj . It has not been theoretically proved, however, that eq. (j2.3j) is the absolute minimum for 
the improved gauge actions. This is checked only numerically 9 . 



2.2 Gauge action 

Our improved action includes plaquette and rectangle loops, and is given by 

Simp [CI = lE W (C,g 2 )2C(C) + \ ]T W X (C, g 2 )2£(C), (2.6) 

with 

C{C) = ReTr[/- U(C)}, (2.7) 

where U{C) is the ordered product of the link variables along loop C contained in the set 
Sq (plaquette) or S\ (rectangular). So and Si consist of all loops of the given shape which can 
be drawn on the cylindrical lattice with volume L 3 x L. The loops involve only the "dynamical 
links" in the sense specified above, and spatial links on the boundaries at xq = and xq = L. 
In particular, rectangles protruding from the boundary of the cylinder are not included. 

One has to choose appropriate boundary weights to achieve the 1-loop level 0(a) improve- 
ment for observables involving derivative with respect to the boundary. Among various choices 
to achieve this, ours is given as follows. 



Wo{C,gl) 



c s (go) for C G P s : Set of plaquettes that lie on one of 

the boundaries, 
cocf (<7q) for C G P t : Set of plaquettes that just touch one 

of the boundaries, 
cq for C G -Pother : otherwise, 
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for C € R s : Set of rectangles that lie completely 

on one of the boundaries, 
Wi(C, <7g) = < ciq^(<? 2 ) for C € i? 2 : Set of rectangles that have exactly 2 (2.£ 

links on a boundary, 
ci for C € -Rother : otherwise, 



with 



cocf( ffQ 2 ) = co(l + cr (1) 5o 2 + 0( 5o 4 )), (2.10) 
dcfCflg) = Cl (3/2 + cf (1) 50 2 + O( 50 4 )), (2.11) 

where coefficients cq and ci of the improved gauge action are normalized such that Co + 8c% = 1. 
We call qf (<?o) and c^(oq) 0(a) boundary improvement coefficients. So far, 1-loop coefficients 
and cjf^ are independent each other. Weight factors Wj(C,g 2 ), which include loop cor- 
rections, becomes the Choice B of |11| in the weak coupling limit. The Choice B achieves the 
tree-level 0(a) improvement and, at the same time, the lattice background field V in ea. ()2.3|) 
satisfies the equation of motion obtained by the variation of dynamical links. 

Incidentally, we discuss the 0(a) boundary counterterm from a different point of view. Here, 
it is assumed that the plaquette loops and rectangle loops which lie completely on the cylinder 
L 3 x L are included in the action, and that each boundary consists of one time slice only. 
As explained in 1 , at order a in the pure gauge theory, there are two possible boundary 
counterterms, a TrjFojifofe} an d a Tr{Fp.iFki}, each of which are summed over the x$ = or 
xq = L hyper plane. Since the latter boundary term vanishes in the case of Abelian constant 
boundary field, in the following, we consider only the former. In this case, we have three 
candidates for 0(a) boundary counterterms which respect lattice symmetries, 

1. spatial sum of time- like plaquette loop that just touches one of the boundaries, 

2. spatial sum of rectangle loop that has exactly 2 links on a boundary, 

3. spatial sum of rectangle loop that has exactly 1 link on a boundary, 

to satisfy one condition, the 0(a) improvement condition. Therefore, for simplicity, we can take 
a trivial weight for the term 3, and we still have one degree of freedom for the choice of the 
boundary terms. At the tree level, however, the background field given in eq, (|2.3|) must satisfy 
the equation of motion 1 , so that one has to take cf = 1 and = 3/2 (Choice B). Since no such 
an extra constraint exists for the 1-loop boundary terms, we can freely set a relation between 
and cjf , which will be given in the next section. 

3 SF coupling and O(a) boundary improvement coefficients 

The SF with the improved gauge action is given by 

Z = e~ r = [ D[U]e- s W, (3.1) 



lr The equation of motion for the plaquette action is trivially satisfied. 
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where S[U] = Si mp [U]. We require the same boundary condition, ea. (|2.2j) . for the link variables 
as in the case of the Wilson plaquette action. In perturbative calculations, there are two main 
concerns to note : one is whether the background field given by eq. ()2.3|) corresponds to the 
absolute minimum of the action, and the other is the gauge fixing. For the latter, we used the 
covariant gauge fixing procedure outlined in [1 _. The former statement is positively proved in 
for Wilson plaquette action. Unfortunately the statement has not been proved yet in the 
case of the improved gauge actions, since the proof in P is not applicable to these cases. In 
however, it has been numerically checked that the background field given in eq. Q2.3|) corresponds 
to the minimum for a large class of improved actions, hence we assume this in our perturbative 
calculations. 

In a neighborhood of the background field V, any link variables U can be parametrized by 

U(x, y) = eyLY>{g aq lx {x)}V(x, /z), (3.2) 
where are quantum fields. The SF coupling is defined through the free energy F in eq. ()3.1j) 



ri=i/=0 



(3-3) 



where V is the derivative with respect to ij. T' Q is a normalization constant 



r' 



90=0 



d_ 

Oi] 



12( - 

a 



CQ 



L 



J2mc)+ E 2£(c) + Cl |E 2£ ( C )+ E 2 ^C) 



Pt 



u=v 



cq sin 27 + sin 7 + 4ci sin 47 + sin 27 



(3-4) 



where 7 is given in appendix TK\ 

P(l) 

Let us discuss the perturbative expansion of the SF coupling. If we require that c t and 
c t satisfy (this is possible by using the last degree of freedom as mentioned in the end of the 
previous section) 

2cf (1) , (3.5) 

(3-6) 



and introduce c t 



(i) 



cocf«+4 Cl cfW 



then one obtains the relation between S[V] and Tq 

1 2a 



S[V] 



1 + L c t 



r + o( 5o 2 ) 



Using the eq. (|3.7|) . the perturbative expansion of g| F (L) is given by 



9sf(L) 



9l 



m 



(i) 



{L/a)gt + 0(gt), 



(3.7) 



(3. 



5 



with 

m£ 1} (L/a) = -^c{ 1} + mS 0) (L/a), (3.9) 

where m^{L/a) is the 1-loop correction to the SF coupling, calculated with the tree-level 0(a) 
boundary coefficients, and the detail of the calculation will be given in the next section. The 
value of the 1-loop coefficient c[ is determined by the condition that the dominant part of the 
scaling violation of m^~\L/a) should be proportional to (a/L) 2 , and then cf ^ and c^ 1 ' are 
uniquely given by ea. (|3.5|) and ea. ()3.6j) . 



$0*0 = ^E eipX % a (P^o), (4.2) 



p 



4 Calculation of the 1-loop coefficient 

In the following, we choose lattice unit (i.e. a = 1). According to the unpublished note ^2j, we 
have used I a (a = 1, 2, • • • , 8) as a basis of Lie algebra of SU(3) in the presence of the background 
field. Their explicit form can be found in ^3]. Decomposing in a basis I a 

a 

the quantum fields are Fourier transformed with respect to spatial momenta as 

1 

%( x ) = ^3E eipXe§(Pfe+ * a(:Co)) 9l(P^o), (4.3) 
p 

where the phase (j) a {xo) is given in appendix El In terms of <7«(a;), the quadratic part of the 
improved gauge action ea. (|2.6|) takes the form 

5 im P = 7l£ £ E^(-P^o)^(p;x ,y )^(p,yo), (4.4) 

P xo, 3/0=0 a 

with the condition 

%(P^o)L 0=0 = 0, for A; = 1,2, 3. (4.5) 
The explicit form of the inverse propagator is given in and also in appendix El 



The 1-loop correction m\ (L) is given by 



m?\L)= ' " 



F'n di] [2 



1 In DetK - In Det A 



(4.6) 

r]=u=0 



where the determinant for the quantum field sector (the first term in the right hand side of the 
equation) is taken with respect to the spatial momentum p, the time xo, the Lie algebra sector 
a and Lorentz index ju. The second term in the right hand of ea. (|4.6|) represents a contribution 
from the ghost sector. Here we will exclusively consider the quantum field sector, since the 
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contribution from the ghost sector to the 1-loop correction is same as in the case of Wilson 
plaquette action. Our boundary condition for the temporal component go is different from that 
in pQ, so that the "non- uniform" contribution in the gauge fixing term remains in the inverse 
propagator K (see appendix |A"|) . 

We evaluated the 1-loop correction rrij (L) numerically for the Iwasaki action (c\ = —0.331, 
c 2 = c 3 = 0) [H], the Ltischer-Weisz (LW) action (ci = -1/12, c 2 = c 3 = 0) |E| and the DBW2 
action {c\ = —1.40686, c 2 = C3 = 0) ^j] i n the range L = 6, • • • ,48. The results are shown 
in Table ^ The computations have been performed by using FORTRAN with the extended 
precision. As a check of our calculation, we have confirmed an independence of the gauge fixing 
parameter and expected symmetries before reducing the amount of the calculation ^7]. We 
have also checked that our code at c\ = reproduces the known result of the Wilson plaquette 
action ^H]. Furthermore, two codes written independently by the two authors have produced 
identical results up to about 30 digits in the range L = 6, • • • , 32. Beyond this range, we have 
used the faster code only. 



5 Analysis and results 

In this section we extract the order 1/L term from the 1-loop correction to^(L) to determine 
the O(a) boundary improvement coefficient c[ . According to the Symanzik's analysis of the 
cutoff dependence of Feynman diagrams on the lattice, one expects that 1-loop coefficient has 
an asymptotic expansion 

00 

m (0) (L) L^oo Y^{A n + B n \nL)/L n . (5.1) 

n=0 

Using the blocking method of |19j . we extracted the first few coefficients Aq, Bq, A%, B\ and 
estimated their errors. 

Some of these coefficients are known or related to other quantities: For example, ^4o' s of 
two different actions are related to the ratio of A-parameters of two actions. If the ultra-violet 
divergences in the SF is removed by the standard renormalization of the coupling constant, 
Bq = 26o, where 60 = ll/(47r) 2 is the 1-loop coefficient of the /3-function in the pure SU(3) 
gauge theory. If the tree- level O(a) improvement is implemented, B\ = must hold. Our main 
result comes from A\: ea. (|3.9|) gives = A\/2. 

We have first verified that our extraction of -Bo an d B\ is consistent with the above expec- 
tation. We have confirmed Bq = 2&o up to 7 digits (Iwasaki), 9 digits(LW) or 4 digits(DBW2), 
while B 1 < 10~ 4 (Iwasaki), < 10" 7 (LW) or < 10~ 2 (DBW2). Since our data give expected values 
of So and B\, we fix Bq = 26o and B\ = by hand, in the blocking procedure to extract Aq 
and A\, whose results for each action are shown in Table |^1 where we have added the result of 
plaquette action [HI El for a later reference. 

As a further check, we extract ^4o' s from the ratio of A-parameters between two schemes X 
and Y, which is given by 

Ay £ 

-± = e «o, (5.2) 



7 



T 
l . 


Iwasaki action 


LW action 


DBW2 action 








0865021015584032 





3843092560841445 


-0 


2542597063902088 


7 





1026697312426737 





4061279685078025 


-0 


2517151449619943 




o 





1171638577366678 





4249279311929165 


-0 


2462340808659547 


Q 

y 





1303628849788211 





4414718008748639 


-0 


2394316692834335 


i n 

1U 





1424404981803593 





4562496675217995 


-0 


2322864496797361 


1 1 
ii 





1535565273026473 





4696045216347286 


-0 


2251762317750280 


1 9 
1Z 





1638489264935022 





4817875167578513 


-0 


2182314272034601 


1 3 
lo 





1734301653189940 





4929883879248706 


-0 


2114993870451222 


1 /I 
14 





1823916642888328 





5033540747374835 


-0 


2049961094379959 


10 





1908085280119190 





5130007428050310 


-0 


1987234947986162 


ID 





1987431354745856 





5220218429952182 


-0 


1926766886112917 


1 7 
1 1 





2062478088988484 





5304936855855431 


-0 


1868475754585161 


lo 





2133668331494581 





5384793988419066 


-0 


1812265016307355 


1 Q 

iy 





2201380038925787 





5460318047790646 


-0 


1758031806343642 


zu 





2265938267004754 





5531955499144914 


-0 


1705672082184823 


91 
Zl 





2327624550322364 





5600087116401930 


-0 


1655083693982648 


99 
zz 





2386684311658198 





5665040280450621 


-0 


1606168213494539 


93 
ZO 





2443332770354413 





5727098525006191 


-0 


1558831969805698 


9/L 
Z4 





2497759696349747 





5786509038428871 


-0 


1512986565129531 


9^ 
ZO 





2550133267993803 





5843488625662739 


-0 


1468549050440217 


9fi 
ZO 





2600603227656853 





5898228494962823 


-0 


1425441882890768 


97 
Z f 





2649303482326239 





5950898137064368 


-0 


1383592748238773 


98 
ZO 





2696354261875951 





6001648495876492 


-0 


1342934304660328 


9Q 

zy 





2741863922040979 





6050614580595203 


-0 


1303403885609065 


3n 





2785930459882241 





6097917633368525 


-0 


1264943186395051 


ol 





2828642794964034 





6143666940320125 


-0 


1227497950252100 


39 
OZ 





2870081858351735 





6187961354134090 


-0 


1191017663623826 


33 





2910321522986861 





6230890581649763 


-0 


1155455266354859 


3/1 





2949429402366960 





6272536278701680 


-0 


1120766879801252 


3^ 
OO 





2987467539279429 





6312972985837565 


-0 


1086911554135926 


3fi 





3024493002264770 





6352268931891717 


-0 


1053851035018207 


37 





3060558404258739 





6390486727199947 


-0 


1021549549113093 


38 
OO 





3095712355291510 





6427683964162490 


-0 


0989973607544213 


3Q 
oy 





3129999859059953 





6463913739632261 


-0 


0959091826148500 


zin 

4U 


n 


316346266152591 1 


n 


64992251 1 1032890 


-0 


0928874761 305482 


41 





3196139558344332 





6533663496047807 


-0 


0899294760096131 


42 





3228066666825220 





6567271024057376 


-0 


0870325823576189 


43 





3259277667231945 





6600086846150990 


-0 


0841943482007244 


44 





3289804017476282 





6632147409439818 


-0 


0814124680962679 


45 





3319675144656546 





6663486700493243 


-0 


0786847677306482 


46 





3348918616375112 





6694136461978401 


-0 


0760091944125418 


47 





3377560294345988 





6724126385966856 


-0 


0733838083775785 


48 





3405624472446620 





6753484286860957 


-0 


0708067748282852 



Table 1: 1-loop coefficient m\ (L) for improved actions. 
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where 

&(p) = 9x(n) + <$x(ti + ---- (5-3) 
A purely numerical number c here is given by 

c = A*-A%, (5.4) 

where Aq or Aq is the expected Aq of the scheme X or Y, respectively. We then find 

"A, 



^imp = A plaq_ 26oln 



'-imp 



Aplaq 



Amp 



(5.5) 



59.05 ± f .0 for Iwasaki action |2T] 

5.29 ± 0.01 for LW action [22] , (5.6) 



Using Ag laq = 0.36828215(13) [HI CHI and the ratio of A-parameters 2 , 



A piaq j 1308 for DB W2 action E 



we obtain the values of A Q mp for each action, which are shown in Table |2] (Aq XP ). We have 
observed the consistency in Aq between previous known results and our calculations. 

With these confidences in our computation, we obtain the main result of our paper, the 
1-loop O(a) boundary improvement coefficient eq. (|3.9|) . which is given by 

c« = cf (1) = Ax/2, (5.7) 

where A\ is also found in Table |2j 

Finally, using our results Aq in Table El we can estimate the ratio of A-parameters between 
the improved action and the plaquette action more accurately, 

Aim P = A imp/A S F = j J_ ^plaq _ ^ 

Aplaq A plaq / A S F I 26 

61.2064(3) for Iwasaki action 

5.292104(5) for LW action . (5.8) 
1273.4(8) for DBW2 action 

6 Conclusion and discussions 

Combining our result of qp for the improved gauge actions and the previous result of for 
the clover quark action |7j, we obtain 



3 f«=cfW/2 = A 1 /2 + n /C f« (6.1) 



for the n f flavors QCD, where cf (1) = 0.0191410(1). 
2 We take the value for DBW2 action from a private note p 
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As a final remark, let us discuss the lattice artifact of the step scaling function (SSF) [2U] 
for various gauge actions. The SSF a(s,u) describes the evolution of a renormalized coupling 
under finite rescaling factor s (say s = 2) 

a(s,u) = g 2 (sL) (6.2) 

u=g 1 (L) 

and it has a perturbative expansion 

a(s,u) = u + 2b ln(s)u 2 + 0(u i ). (6.3) 

This SSF a(s,u) in the continuum theory is obtained by the continuum limit of the lattice SSF 
£(s, u, l/L): 

cr(s,u) = Jim T,(s,u,l/L). (6-4) 

Therefore we can estimate the lattice artifact of the SSF in our perturbative calculation. We 
define the relative deviation 5(s,u, l/L) and expand it as 

6(2, u, l/L) = n^l/L)-a{2,u) = ^ + 
a(2,u) 

(k) 

where we have set that s = 2 and 5{ (2, l/L) is the 1-loop coefficient. Here k denotes the 
degree of the improvement for the boundary coefficient: the tree (1-loop) value is used for k = 
(k = l). 

The manifest form of 5^(2, l/L) is given by 

5{ k) (2,l/L) =m[ k) (2L) -m[ k) (L) -26 ln(2), (6.6) 

and the results, including data of the plaquette action [HIE] for comparison 3 , are given in Table 
121 for each gauge action. Figure ^ (Figure I^J) shows that the 1-loop deviations with the tree-level 
(1-loop level) 0(a) improved boundary coefficient vanish roughly linearly(quadratically) in l/L 
as expected. As is evident from Figure 1 and Figure 2, at 1-loop level, the lattice artifact for the 
renormalization group improved action (Iwasaki or DBW2) is comparable to or larger than that 
for the plaquette action, while the LW action is the least affected by the lattice cutoff. However, 
one can not conclude that the LW action is the best choice for numerical simulations, where the 
lattice artifacts of higher orders in u or a may not be negligible. 
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We have added data for the plaquette action in the range L — 17, • ■ ■ , 24 
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plaquette action 


Iwasaki action 


LW action 


DBW2 action 


^exp 
k 


0.36828215(13) 
-0.17800(10) 


-0.2049015(4) 
-0.1999(24) 
0.30360(26) 


0.136150567(6) 
0.13621(26) 
-0.005940(2) 


-0.62776(8) 
-0.6159 
0.896(45) 



Table 2: The coefficients of asymptotic expansion Aq, A± for various gauge actions. The values 
for plaquette action are taken from [SJ 





plaquette action 


Iwasaki action 


LW action 


DBW2 action 


L 


4°) 
















6 


0.01089 


-0.00394 


-0.01922 


0.00608 


0.000911 


0.000417 


-0.061 


0.014 


7 


0.01004 


-0.00268 


-0.01684 


0.00484 


0.000659 


0.000236 


-0.050 


0.014 


8 


0.00918 


-0.00194 


-0.01499 


0.00399 


0.000527 


0.000156 


-0.043 


0.013 


9 


0.00841 


-0.00148 


-0.01356 


0.00330 


0.000441 


0.000111 


-0.038 


0.011 


10 


0.00773 


-0.00117 


-0.01241 


0.00277 


0.000379 


0.000082 


-0.035 


0.010 


11 


0.00714 


-0.00095 


-0.01146 


0.00235 


0.000333 


0.000063 


-0.032 


0.009 


12 


0.00663 


-0.00079 


-0.01064 


0.00201 


0.000296 


0.000049 


-0.030 


0.008 


13 


0.00618 


-0.00066 


-0.00994 


0.00174 


0.000268 


0.000039 


-0.028 


0.007 


14 


0.00579 


-0.00057 


-0.00932 


0.00152 


0.000244 


0.000032 


-0.026 


0.006 


15 


0.00544 


-0.00049 


-0.00878 


0.00134 


0.000224 


0.000026 


-0.024 


0.006 


16 


0.00513 


-0.00043 


-0.00830 


0.00119 


0.000207 


0.000022 


-0.023 


0.005 


17 


0.00486 


-0.00038 


-0.00787 


0.00106 


0.000193 


0.000018 


-0.022 


0.005 


18 


0.00461 


-0.00034 


-0.00748 


0.00095 


0.000181 


0.000016 


-0.021 


0.004 


19 


0.00438 


-0.00030 


-0.00713 


0.00086 


0.000170 


0.000013 


-0.020 


0.004 


20 


0.00418 


-0.00027 


-0.00681 


0.00078 


0.000160 


0.000012 


-0.019 


0.004 


21 


0.00399 


-0.00025 


-0.00652 


0.00071 


0.000151 


0.000010 


-0.018 


0.003 


22 


0.00382 


-0.00022 


-0.00625 


0.00065 


0.000144 


0.000009 


-0.017 


0.003 


23 


0.00367 


-0.00020 


-0.00601 


0.00059 


0.000137 


0.000008 


-0.017 


0.003 


24 


0.00352 


-0.00019 


-0.00578 


0.00054 


0.000131 


0.000007 


-0.016 


0.003 



Table 3: The deviations for various gauge actions. 



11 



0.02 



0.01 


-0.01 
-0.02 
-0.03 
-0.04 
-0.05 
-0.06 
-0.07 



□ □ 



□ □ a 



□ □ 



AAA A A A A A A A A A & £ A_ 



( ^oo 0o 



O g 



plaquette action □ 

Iwasaki action o 

LW action & 

DBW2 action o 



_i i_ 



0.02 0.04 0.06 0.08 0.1 

1/L 
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Figure 1: The relative deviations of the lattice SSF from the continuum one at 1-loop for 



various gauge actions with the tree-level O(a) improved boundary term. One can find that S\ 
for various gauge actions vanish roughly linearly in 1/L. 
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Figure 2: The same quantities as in Fig. ^ but with the 1-loop level O(a) improved boundary 
term. One can find that 6-j/ for various gauge actions vanish roughly quadratically in 1/L. 
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A Inverse propagator 



Here, we give the explicit form of the inverse propagator. We choose the lattice unit (i.e. a = 1) 
and set T = L. 



Plaquette 



K = c K (0) + ciK« + A if (gf) . 



^OO^P; X 0, 2/0 ) = R a Sx ,y S a (P, X ) ■ S a (p, X + 1), 



^!o )a (P; x o,2/o) 



Klf a {p;x ,y ) 



Rectangle 

Wdbc = 3/2 for the choice B |TTj . 



iR a 



S x ,yo s k(.P' x O + 1) - ^ -l,j/o4(P>2/o) 



- K ko a (p;yoiX ), 



J %o,yo 



5 M s a (p,x ) 2 - st(p,x )s?(p,x ) 
2C 5 Xo>yo — R {S Xo +i,y + S XQ — 1,2/ 



^oo )a (P; x o,2/o)ofcfe 



J xo,yo 



1 + (w^dbc - 1)(^ ,0 + ^ ,L-l) 



x4i^^sin(> a (xo) + p m ) sin((j) a (x + 1) +p m ), 



^oo )a (p; 2; o,2/o)oofc 



(1 - tf:eo,I.-l) s m(P> x o)Sm(P> x + 2) 



+ (1 - *x ,o)flm(P» x - l>m(P, ^0 + 1) 
-l,2/0 S m(P> x o 
,yo s m 

(p,y - l)s^(p,y + 1) 



Kj^ a (p;xo,yo)okk 



iR a 2 2c a k (p,x ) 



(A.l) 

(A.2) 

(A.3) 
(A.4) 



(A.5) 



(A.6) 



(A.7) 



+ ( W dbc ~ l)S XQtL -l}sm((j) a (xQ + 1) +p k ) 



5x -l,y { 1 + ( w dbc ~ l)4o,l}sin(0 a (2/o) + Pfc) 



, (A.8) 



K^ a (p;x ,yo)ook = i^{s%(p,xo + 2)[(l-S Xo ,L-i)8x ,yo + s xo+i 



,'/oJ 
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- st(p,x - 2)[(1 - 5 xoA )5 xo _ hyo + <y X o-2,vo]} > ( A - 9 ) 

(A.10) 



-^feo )a (p;yo,ajo), 



-^S )a (p;a;o,yo) 



Ofcfc 



<5jwcjfc(p,a;oK(p,j/o) 

+ (w d6c - l)fe^ ,j /0 Cfc(p,X ) 



IC^xq^q ^2(^x0+1,2/0 + ^x -l,2/o, 



4 ,i(C 2 a 4(p,x )-i5f4(p,xo)) 



^fc/ )a (P;^0,yo)other S 



XQ.J/O 



+ <W-i(C 2 a c£(p, x ) + iS$s%(p, so)) 
*fci S S ™(P' x o) 2 (cfc(p, x ) 2 + C(p, x ) 2 ) 

m 

s£(p, x )s?(p, x )(cg(p, x ) 2 + cf (p, x ) 2 ) 



,(A.ll) 



(2 - 5 XOj l - 5a;o,L-l)C2^o,J/0 ~~ ^2(^x0+2,2/0 + $r -2. 



;vo / 



.(A.12) 



Gauge fixing term 



^oo f) °(p;a?o,W)) 



^o f)a (p;^o,yo) 



^f )a (p;^o,yo) 



4f )a (p;xo,yo) 



— 2(5 XOiJ/0 ^xo+i,ko ^xo-1,2/0 



-^o,2/o[^o,o(l - Xa5p,o) + *x ,L-l] 
= -Wfc(P,iCo)[*xo,VO -*x -l,»d]. 



-^■o f)a (p;yo,^o), 



Coefficients 



4(p,x ) = 2sin[(p fc + a (x o ))/2], 
ct(jP,xo) = 2cos[(p fc + a (x o ))/2], 

4>a (xq) = -<t>a(xo), 

s~ia s~id rya r>a rya 

u — o , o — — O , JX — ft , 

Ca /-fa c*a na pa pa 

2 — 2 , J'2 — 1 n 2 — Il 2 i 

Xa = Xa= (0,0,1,0,0,0,0,1), 



(A.13) 
(A.14) 
(A.15) 
(A.16) 

(A.17) 
(A.18) 
(A.19) 
(A.20) 
(A.21) 
(A.22) 



15 



where 1 = 2, 4 = 5, 6 = 7, and vice verse. For diagonal part, 3 = 3, 8 = 8. 
Lists of coefficients 



a 




S a 


1,4 


| (cos 27 + cos 7) 


— i^(sin 27 + sin 7) 


3,6 


cos 7 





8 


5 (2 cos 27 + cos 7) 






Table 4: C a and S a for SU(3). 



a 


<t>a{xo) 


i? a 


1 




cos ^ 


4 


-37x0 + K»/[§ + 


COS ^ 


3 





cos 7 


6 


r(2^-i) 


cos 7 


8 





|(2cos 27 + cos 7) 



Table 5: 4> a (x ) and R a for SU(3). 



a 


it 2 




qa 


1,4 


cos 7 


i?2 COS 37 


—iR\ sin 37 


3,6 


cos 27 


R a 2 





8 


i(2cos 47 + cos 27) 








Table 6: R%, C$ and S£ for SU(3). 
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